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ABSTRACT 

We continue the analysis of non-Gaussianities in the CMB by means of the scahng 
index method (SIM, Rath, Schuecker & Banday 2007) by applying this method on 
the single Q-, V-, W-bands and the co- added VW-band of the 5-year data of the 
Wilkinson Microwave Anisotropy Probe (WMAP). We compare each of the results 
with 1000 Monte Carlo simulations mimicing the Gaussian properties of the best fit 
ACI^M-model. Based on the scaling indices, scale-dependent empirical probability 
distributions, moments of these distributions and x^-combinations of them are cal- 
culated, obtaining similar results as in the former analysis of the 3-year data: We 
derive evidence for non-Gaussianity with a probability of up to 97.3% for the mean 
when regarding the KQ75-masked full sky and summing up over all considered length 
scales by means of a diagonal x^-statistics. Looking at only the northern or southern 
hemisphere of the galactic coordinate system, we obtain up to 98.5% or 96.6%, respec- 
tively. For the standard deviation, these results appear as 95.6% for the full sky (99.7% 
north, 89.4% south) and for a x^-combination of both measurements as 97.4% (99.1% 
north, 95.5% south). We obtain larger deviations from Gaussianity when looking at 
seperate scale lengths. By performing an analysis of rotated hemispheres, we detect 
an obvious asymmetry in the data. In addition to these investigations, we present a 
method of filling the mask with Gaussian noise to eliminate boundary effects caused 
by the mask. With the help of this technique, we identify several local features on the 
map, of which the most significant one turns out to be the well-known cold spot. When 
excluding all these spots from the analysis, the deviation from Gaussianity increases, 
which shows that the discovered local anomalies are not the reason of the global de- 
tection of non-Gaussianity, but actually were damping the deviations on average. Our 
analyses per band and per year suggest, however, that it is very unlikely that the 
detected anomalies are due to foreground effects. 

Key words: cosmic microwave background - cosmology: observations - methods: 
data analysis 



1 INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite, launched in June 2001, measures the temperature 
anisotropy of the cosmic microwave background (CMB) ra- 
diation with surpassing accuracy, hence providing the best 
insight on the beginnings of our universe until now. From 
the first data release on, many investigations were made 
concerning the Gaussianity of the CMB, since such analyses 
give information about the nature of the primordial den- 
sity fluctuations, which are the seeds of those temperature 
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anisotropics. The statistical properties of the density fluctu- 
ations are again an important observable for testing cosmo- 
logical models, especially models of inflation. Standard infla- 
tionary models predict the temperature fluctuations of the 
CMB to be a Gaussian random field which is isotropic and 
homogenous ( |Guth|198ll |Linde|1982| [Albrecht Steinhard 
1982). Still, there also exist more complex models that allow 
non-Gaussianity in a s cale-independend (Linde Mukhanov 
19971 |Peebles| [19971 [Bernardeau Uzan||2002| |Bartol" 



Mata rrese Riotto |2002| |Lyth W ands "20Q2| |Lyth, Un- 



garelli &: Wands [2003 Acquaviva et al. 2003 ) or in a scale- 



dependend way ([Garriga Mukhanov ||1999[ Armendariz- 
Picon, Damour Mukhanov [1999 ). For a detailed overview 
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of the different models and a more specific survey on scale- 



I 



dependend ones, see Bartolo et al. ( 2004 ) and Lo Verde et al 



(2008), respectively, as well as enclosed references. In addi- 
tion, topological defects like cosmic strings can induce local 
non-Gaussianities and influence the power spectrum (|Kaiser| 



&: Stebbins [T oMI'Bo uchet, Bennett, &; Stebb ins'1988^'Turok' 
Spergel 119 90 i Turo k |1996||Jeanneroti^ cher & Sakellar- 



iadou ||2003| . Considering this plethora of possible physical 



mechanisms, which may induce non-Gaussianity, studies of 
Gaussianity of the CMB are strongly required for testing 
predictions of fundamental physical theories. By compar- 
ing the results with theoretical predictions, we can evaluate 
which model e.g. of inflation can be accepted or rejected. 

Non-Gaussianity implies the presence of any higher or- 
der correlations. Therefore, a concrete description of the 
characteristics of non-Gaussianity is not possible, and one 
has to state that it can occur in various forms. One can carry 
out a global analysis to search for deviations from Gaussian- 



ity dKomatsu et al. |2003 


2009||Chiang et al. |2003||Chiang, 


Naselsky & Coles ||2007| 


Coles et al. ||2004| |0 'Dwyer et al. 


112004' 'Eriksen et al.||2007| 


Smith, Senatore & Zaldarriaga 


2009; Rudjord et al. ||2009 


). But one can also concentrate 



on more specific investigations (this is most often performed 
in addition to a general analysis) , whereas we want to point 
out the following two: 

Investigations concerning asymmetries in the CMB 
data were accomplished with linear (Eriksen et al. | 2004a, 



Bonaldi et al. |2007||Bernui 


|2008||Bernui & Hipolito-Ricaldi| 


|2008 


Hoftuft et al. ||2009 


) as well as non-linear methods 



Hansen et al. ||2004a| [Land &: Magueijo ||2005| |Pietrobon| 
et al. |2008||Rath et al. |2009| ). With those methods, studies 



of the differences between the northern and southern hemi- 
sphere of the galactic coordinate system, naturally given by 
the absent region of the out masked galactic plane, as well 
as a search for a preferred direction of maximum asymme- 
try were performed. In almost all investigations, significant 
asymmetries between the north and the south were detected. 
Thereby, it depended on the type of analysis, which hemi- 
sphere featured the larger deviations from Gaussianity and 
which hemisphere agreed better with the standard model. 
The preferred direction of maximal disparity was in most 
investigations found to lie close to the ecliptic axis. 

Local features are another particular form of non- 
Gaussianity being of growing importance, e.g. for the search 
of topological defects like cosmic strings. Since the first de- 
tection of the famous cold spot by|Vielva et al. (2004), many 



investigations tried to find new, or re-detect known spots 



with various methods (Mukherjee & Wang^ 2004 Cayon, 


Jin & Treaster |2005| 


Cruz et al. |2005||2006||McEwen et al. 


20051 12006| 20081 |Vie 


Iva et al. |2007| |Pietrobon et al. |2008| 


Gurzadyan & Kocharyan |2008| |Gurzadyan et al. |2009|. In 



doing so, several significant spots have been detected up to 
now. 



In Rath et al. (2007), all mentioned investigations were 
accomplished by applying, for the first time, the scaling 
index method on the WMAP 3-year data. In this paper, 
we continue these analyses by applying the scaling index 
method on the WMAP 5-year observations. We search for 
global non-Gaussianities and asymmetries in the data and 
use a modified approach to detect local features. 



This paper is structured as follows: In Section 2 we 
present the preprocessing of the WMAP data and the 
modality of creating the simulations. In Section 3, the scal- 
ing index method is introduced as well as a technique to cope 
with boundary effects. With these requisites, we are ready 
to perform our calculations, whose results are presented in 
Section 4. In this chapter, we first discuss the global inves- 
tigations as well as asymmetries and focus on local features 
later on. All these findings are summarised in Section 5. Fi- 
nally, we draw our conclusions in Section 6. 



2 WMAP DATA AND SIMULATIONS 

For our investigations we use the Q-, V- and W-band 
five-year-data of the WMAP-satellite as it is provided by 
the WMAP-Tean|5 We work with the foreground-reduced 
maps, which use the Foreground Template Model proposed 
in iHinsh aw et al. | ( [2007 ) and Page et al. | (|2007|) for fore- 
ground reduction. To obtain a co-added VW-map as well as 
single V-, W- and Q-maps, we accumulate the differencing 
assemblies Ql, Q2, VI, V2, Wl, W 2, W3, W4 via a noise- 
weighted sum |Bennett et al. ||2003[ ): 



no, 



(1) 



In this equation, A characterises the set of required 
assemblies, e.g. for the co- added VW-map A = 
{V1,V2,W1,W2,W3,W4}. The parameters 6 and cor- 
respond to the co-latitude and the longitude on the sphere, 
while the five-year noise per observation of the different as- 



semblies, given by Hinshaw et al. (2009), is denoted by ctq 



We decrease the resolution of the maps to 786432 pixels, 
which equ als to Ngjde — 256 in the employed HEALPix- 
softwara^ (Gorski et al. 2005) and cut out the heavily 



foreground- affected parts of the sky using the KQ75-mask 
(Gold et al. [2009 ), which has to be downgraded as well. We 
choose a conservative downgrading of the mask by taking 
only all pixels at Nside — 256 that do completely consist of 
non- mask-pixels at Nside = 512. All downgraded pixels at 
Nside — 256, for which one or more pixels at Nside — 512 
belonged to the KQ75-mask, are considered to be part of the 
downgraded mask as well. In doing so, 28.4% of the sky is 
removed (see upper left part of figure [l]). Finally, we remove 
the residual monopol and dipol by means of the appropriate 
HEALPix routine applied to the unmasked pixels only. 

To accomplish a test of non-Gaussianity, we also need 
simulations of Gaussian random fields. We create 1000 sim- 
ulations for every band and proceed hereby as follows: We 
take the best fit AC DM power spectrum Ci, derived from 
the WMAP 5-year data only, and the respective window 
function for each differencing assembly (Q1-Q2, VI- V2, Wl- 
W4), as again made available on the LAMBDA- website^ . 
With these requisites, we can create Gaussian random fields 
mimicing the Gaussian properties of the best fit ACDM- 
model and including the WMAP-specific beam properties 
by convolving the Ci 's with the window function. For every 
assembly, we add Gaussian noise to these maps with a par- 
ticular variance for every pixel of the sphere. This variance 



^ http://lambda.gsfc.nasa.gov 
^ http://healpix.jpl.nasa.gov 
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Figure 1. The two plots on the left hand side illustrate the original 5-year WMAP-map of the co-added VW-band (above) and the 
related colour-coded a-response (below). The equivalent plots for the mask-filling technique are arranged on the right hand side. These 
maps (and all following ones) are shown in a conventional scheme, namely the Mollweide projection in the Galactic reference frame with 
the Galactic Centre at the centre of the image and the longitude increasing from there to the left-hand side. 



depends on the number of observations Ni{0,(l)) in the re- 
spective direction and the noise dispersion per observation, 
(70,2. After this procedure, we summarize the Q-, V- and 
W-bands and the co-added VW-band using equation ([T]), 
decrease the resolution to Nside = 256, cut out the KQ75- 
mask and remove the residual monopol and dipol, just as 
we did with the WMAP-data. 



Radius 


0.025 


0.050 


0.075 


0.100 


0.125 


lhj2] 


[83,387] 


[41,193] 


[28,129] 


[21,97] 


[17,77] 


Radius 


0.150 


0.175 


0.200 


0.225 


0.250 


[hj2] 


[14,65] 


[12,55] 


[10,48] 


[9,43] 


[8,39] 



Table 1. The angular scales corresponding to the position of the 
90% (^i) and the 10% (£2) weighting in the scaling index formula 
when using a given scale parameter r. 



3 WEIGHTED SCALING INDEX METHOD 
3.1 Formalism 

We perform our investigations using the scaling index 
method (SIM) (Rath et a l. |2002||Rath k Schuecker |2003| ), 
which enables a characterisation of the structure of a given 
data set. It has already been used in time series analysis of 
active galactic nuclei (AGN) ( Gliozzi et al. [ [2002 



Gliozzi, 



Papadakis & Rath 2006) as well as in structure analysis 
for 2D and 3D image data, e.g. in Jamitzky et al. (2001); 



Monetti et al. (2003); Rath et al. (2008). In the following, 



we only present a short overview of the calculation of scaling 
indices. For a more detailed formalism of using the SIM in 



CMB analysis we refer to Rath et al. (2007). 



The fluctuations of the temperature maps are charac- 
terized by the values of the pixelised sky of a sphere with 
radius R. We transform this representation to variations in 
the radial direction around the sphere by applying a jitter 
depending on the intensity of the fluctuation. Thereby we 
obtain a point-distribution in the three-dimensional space. 
Thus, given Npix as the number of pixels on the sphere, the 
value of every pixel (^i, (/)i), i = 1, Npix corresponds to a 
vector Pi in the three-dimensional space. We then define for 



every point pi its scaling index by 



a{p^ 



(2) 



ENpix 



where dij denotes the euclidian distance measure 



dij = \\pi -Pjh 

between the points pi and pj, while r characterizes a scale 
parameter. This parameter does not draw a clear-cut line 
between the pixels that are included in the calculations and 
those that are excluded; it rather influences how each sin- 
gle pixel is taken into consideration for the calculation, in 
relation to its distance from the center pixel: For lower r, 
only the closest pixels are important in the calculation of 
a{pi,r), whereas for larger r, the farther distant pixels are 
considered as well, even though with a still lower weight 
than the close pixels. In our study, we use the ten scales 
n = 0.025, 0.05, 0.25, i = 1, 2, 10 and the radius R = 2 
for the sphere. Table ^ shows for each radius r the corre- 
sponding angular scales at the position of the 90% and the 
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10% weighting, thus giving an estimate on how the r- values 
relate to ^-bands in Fourier space. 

The value of a characterises the structural components 
of a point distribution. For example, points in a cluster-like, 
filamentary or sheet-like structure lead to a ~ 0, a ~ 1 or 
a 2, respectively. A uniform distribution of points results 
in a 3, while points in under dense regions in the vicinity 
of point-like structures, filaments or walls have a > 3. 

On the basis of these scale-dependent a- values, we com- 
pute simple measures such as moments and empirical prob- 
ability distributions. We make use of the following scale- 
dependend statistics, namely the mean, the standard devi- 
ation and a diagonal x^-statistics, to compare the results of 
the original WMAP data with the results of the simulations: 



2 



1/2 



2 



Mi{rk) - {Mi(rk)) 



(3) 
(4) 
(5) 



where Mi(rfc) = (a(rfc)), M2(rfc) = (^a(rk) ^ denotes 
the number of pixels in consideration. For all analyses we 
will only consider the non-masked pixels of the full sky or of 
(rotated) hemispheres, as it will be outlined in Section [ 



ionJO 
(2004b| 



Note that we follow the reasoning of Eriksen et al 
and choose a diagonal and not the full x -statistics involv- 
ing the inverted cross-correlation matrix, because also in 
our case the moments are highly correlated leading to high 
values in the off-diagonal elements of the cross-correlation 
matrix. Therefore, the matrix would converge very slowly 
and numerical stability would not be given. If however the 
chosen model is a proper description of the data, any combi- 
nation of measures should yield statistically the same values 
for the observations and the simulations. 

To obtain scale-independent variables as well, we also 
use three diagonal x^-statistics, derived from {a) and aa, 
which sum over all utilised lenght scales r: 



2 

X(a) 



2 

X(q;) ,o 



E 

k = l 

2 Nr 

EE 



Mi(rfc)-(Mi(rfc)) 



M2{rk) - {M2{rk)) 



M.(rfc)-(M.(rfc)) 



(6) 

(7) 
(8) 



Hereby denotes Mi(rfc) = (^(^fe)) and M2{rk) = cra(rfc)- 
The number of different scale parameters r is named Nr. 
Throughout all subsequent investigations, Nr equals ten. 

Finally, to be able to access the degree of difference be- 
tween the data and the simulations and hence a degree of the 
non-Gaussianity of the data, we use the cr-normalised devia- 
tion of the WMAP results of the above-mentioned statistics: 



M - (M) 



(9) 



where in this case M refers to one of the variables defined 
in equation ([3| to (|8| respectively. M itself is calculated by 



Figure 2. A slice of the three-dimensional representation of 
the VW-band WMAP data, illustrated as a XjZ-projection of all 
points with \y\ < 0.05. The upper plot illustrates the original, 
the lower the mask-filling method. The black circles indicate the 
scaling range r = 0.2. 



using the WMAP data, while its moments result from the 
simulations. Note that we pass on the absolute value in this 
general definition to obtain positive as well as negative devi- 
ation. Although we will use the absolute value in the global 
investigations, the sectioning into positive and negative devi- 
ation is useful for the analysis of north-south asymmetry by 
means of rotated hemispheres in chapter [4.1 1 It also allows 
a better interpretation of the character of difference, since 
e.g. a higher mean of the scaling indices implies a more 'un- 
structured' arrangement of the 'pixel cloud' and vice versa. 
Similarily, a higher standard deviation of the indices indi- 
cates a larger structural variability. 

In the tables, we also included the fraction p of the 
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simulations that have higher (lower) values than the data in 
terms of the respective calculated statistics. This percentage 
corresponds to a empirical significance level of the null hy- 
pothesis that the observation belongs to a Gaussian Monte 
Carlo ensemble. 



3.2 Coping with boundary effects 

The regions in the direction of the galactic plane as well as 
many small spots all over the WMAP map are masked out 
since they represent heavily foreground- affected areas which 
would not allow a reasonable analysis of the intrinsic back- 
ground fluctuations. But this operation spoils the results of 
the scaling index method: Instead of a more or less uniform 
distribution, the a- values in the regions around the mask 
now detect a sharp boundary with no points in the masked 
area, into which the scaling regions extends (see figure |2| . 
This results in lower values of a. The effect can clearly be 
seen in the a-response of the masked VW-band WMAP-data 
in the lower left corner of figure [l] A solution to this problem 
is to fill the masked areas with suitable values, that prevent 
the low outcome at the edges of the mask. We accomplish 
this by filling in (nearly) white Gaussian noise with adjusted 
parameters. This is performed by applying the following two 
steps: 

At first, we fill the masked regions with Gaussian noise, 
whose standard deviation for each pixel corresponds to the 
pixel noise made available on the LAMBDA-websitej5 

Here, (J{e,ct)) denotes the pixel noise of the pixel which is lo- 
cated in the direction (^, 0). Then, we scale the expectation 
value and the variance as a whole to the empirical mean firem 
and variance cr^em of the remaining regions of the original 
temperature map: 



Full Sky 



.fe(^,0) 



2^"^ Tmaski^i 0) + Mr 



with 



f^rem 



2 

^raask 



Nm 



where 71 and A4 stand for the non-masked and masked re- 
gion of the map respectively, and N-jz as well as Nm denote 
their number of pixels. Thus, we filled the mask with (nearly) 
white Gaussian noise whose mean and standard deviation 
equal the respective terms of the remaining map, whereby 
the spatial noise patterns are preserved. 

With this filling technique, we obtain a complemented 
data set instead of just excluding the masked regions. Figure 
|2]shows a slice of the three-dimensional representation of the 
temperature fluctuations of both techniques, representing a 
2D projection of the 3D point distribution used for the cal- 
culation of the scaling indices. The center region of the filled 

^ http://lambda.gsfc.nasa.gov 




Figure 3. The probability distributions P(a) of the scaling in- 
dices for the WMAP data (dark lines) and for 50 simulations 
(fainter lines) by using the scale parameter r = 0.2, computed 
for the original (red) and the mask filling method (blue). The 
upper histogram shows the distribution of the full sky data set, 
while the middle and the lower ones show the distribution of the 
northern and southern hemisphere respectively. 



sphere now highly resembles the appearence of the remain- 
ing area, although a more uniform arrangement is visible. 
In the form of a moUweide projection, the filling method as 
well as the corresponding a-response are displayed in the 
right column of figure [l] 

The filling strategy shows obvious success in the adjust- 
ment of the scaling indices map (see the lower right panel 
of figure [T]): The white noise leads to higher values in the a- 
response for the pixels close to the mask as compared to the 
masking method. The regions around the edges of the mask 
feature now a-responses that match far better the values 
of the remaining regions. Still these a- values are calculated 
with the help of an artificial environment, but now the con- 
tortions are lower compared to the original approach. Since 
we apply this method to both the original WMAP data as 
well as to the simulations, the now smaller systematic errors 
in the a-calculation for the 'edge'-pixels are the same for 
both kinds of maps. Thus any significant deviation found in 
the WMAP data is due to intrinsic effects. 

The most important advantage of the filling strategy 
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Figure 4. The histograms of the mean values {a(rk)) of the 
1000 simulations for the five different scale parameters r2, r4, tq, 
rg and rio, calculated for the full sky. The red lines denote the 
corresponding results of the WMAP data, while the black and 
grey lines characterise the average over all simulations and the 
Icr regions, respectively. 



arises if one considers local features in the CMB map: For 
a search of local anomalies (e.g. cold spots), the filled map 
provides a better underlying than the original map: In a 
point distribution, spots as well as points at the border of the 
distribution show similarly low a- values. Thus, by cutting 
out the mask, it is difficult to decide whether a detected local 
feature really exists or whether it originates from a masked 
area nearby. But by using the filling technique, there is no 
longer an edge between the masked and the non-masked 
regions, and anomalies in the Gaussian noise of the mask 
are highly improbable. Thus, any detected local feature must 
then originate from the data itself. Considering the amount 
of masked areas outside the galactic plane, this technique 
describes a considerable improvement for the whole sky. Due 
to these advantages we will only use the mask-filled maps in 
chapter |4.2| 



4 RESULTS 

4.1 Band-wise and co-added map analysis 

In figure [s] the empirical probability densities P{a) of the 
scaling indices (calculated with r = 0.2) are displayed for 
the WMAP data and for the simulations, evaluated for the 
original and the filling method from chapter [3.2[ For clarity 
reasons we only used the first 50 simulations in these plots. 
For both methods, a shift of the WMAP data to higher val- 
ues can be detected, that becomes particularly apparent in 
the northern hemishpere of the galactic coordinate system. 
This indicates a more 'unstructured' arrangement of the un- 
derlying temperature fluctuations of the CMB data in com- 
parison to the simulations. In addition, the histograms of 
the simulations are slightly broader and therewith contain- 
ing a larger structural variability than the one of the WMAP 
data. 

Comparing the non-filling and the filling method, the 
histograms of the latter feature a higher maximum as well 
as higher values for large a, but lower probabilities for 
a G [2.0,2.5]. The obvious reason for this shift is the fact 
that the filled mask does not reduce the a-values of its sur- 
roundings as it was the case with the former method. Now, 
the outcome of these regions is influenced by the white noise 
and is therefore allocated at higher values. 

If we focus on the mean values {a{rk)) of the scaling 
indices, and compare the results of the simulations with the 
WMAP data, the above mentioned shift to higher values 
becomes yet clearer. This can be seen in Figure ^ where 
the distribution of {a(rk)) for the simulations as well as the 
data is displayed for the flve different scale parameters r2, 
r4, re, rg and rio- These results were obtained using the full 
sky. For all applied scales, the distance between the average 
over all simulations and the result of the original WMAP 
map is notably similar. If we perform this analysis for the 
northern hemisphere only, the deviations of the original data 
as compared to the simulations become significantly larger. 

Both the shift to higher values of the WMAP data in 
comparison to the simulations as well as its broader density 
are reflected in the cr-normalised deviations S{r) of the scale- 
dependent statistics of the equations ([s]) - (jsj : The former is 
reflected in the mean and the latter in the standard deviation 
(and therefore both aspects in the diagonal x^-statistics). 
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Figure 5. The cr-normalised deviations S{r) of the statistics of 
the equations (|3| - ([5| in absolute values for the VW-band, plot- 
ted as a function of the scale parameter r. The lines with "+" 
denote the mean, "*" the standard deviation and the boxes the 
-combination, each for the original (red) and the mask filling 
method (blue). As in figure [s] the upper diagram shows the re- 
sults of a full-sky analysis, while the middle and the lower ones 
show the results when only concerning the northern or southern 
hemisphere respectively. 



Figure [5] shows these deviations for the coadded VW-band 
as a function of the scale parameter r for the original and 
the mask filling method, while figure |6] displays only the 
latter method, but for the three single bands. As above, the 
results are illustrated for the full sky as well as for the seper- 
ate hemispheres. The shift to higher values of the WMAP 
data in the northern hemisphere in figure [S] appears now as 
an increased 5'(r), especially for higher scales (r ^ 0.125), 
where the deviations of the two moments range between 2a 
and 3.5(7, and the x^-combination nearly reaches a 6cr-level. 
In the southern hemisphere, only the lowest scales show a 
namable S{r). On larger scales, no signatures for deviations 
from Gaussianity are identified. Looking at the single bands 
Q, V and W, the overall qualitative behaviour of the images 
is quite similar, while the cr-normalised deviations itself are 
slightly lower in most cases. A remarkable fact is the ap- 
pearance of the highest S{r) (5.2cr for the x^-combination in 
the northern hemisphere at r = 0.15) in the frequency band 
which is considered to be the least foreground-contaminated 
one, namely the V-band. Comparing the co- added VW-band 



Figure 6. Same as figurejsjbut applied to the Q- (red), V- (green) 
and W-band (blue). Only the results of the mask-filling method 
are shown while the original method is left out. 



of the original approach and of the mask filling method, 
the cr-normalised deviations of the mean are almost iden- 
tical. The standard deviation of the latter method in com- 
parison to the former one shows a slightly lower S{r) for 
higher scales, which is also reflected in the graph of the 
combination, yet the profile remains the same. 

We also calculated the cr-normalised deviations S and 
the percentages p of the simulations with higher (lower) re- 
sults of the scale-independent diagonal x^-statistics from the 
equations ([6| to ([8| , which are listed in table [2] Although 
the results are damped by a few unimportant scales, high 
deviations are still found, particularly in the northern hemi- 
sphere. For a better comparison to seperate scale lengths, 
the respective results of the scale-dependent statistics Q to 
([5]) are listed in table [s] for which we used the single scale 
r = 0.2. 

In general, all occuring characteristics of the figures |3] 
and [5] match the findings of the analysis of the WMAP 3- 



year data in Rath et al. (2007). This indicates that the 



results are not based on some time-dependent effects. Since 
the 5-year data features lower error bars than the 3-year 
data, it is also improbable that both results are induced by 
noise effects only. 

Evidence for north-south asymmetry in the WMAP 
data was already detected using the angular power spec- 
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Full Sky 


Northern Sky 


Southern Sky 




Full Sky 


Northern Sky 


Southern Sky 




{S/%) 


{S/%) 


{S/%) 


(a{0.2)): 


{S/%o) 


{S/%o) 


{S/%o) 


VW (original) 
VW (mask-filled) 
Q (mask-filled) 
V (mask-filled) 
W (mask-filled) 


2.2 / 97.3 
2.2 / 97.3 
2.1 / 97.4 

2.0 / 97.4 

2.1 / 97.5 


2.7 / 98.5 
2.7 / 98.3 
2.7 / 98.3 
2.6 / 98.1 
2.6 / 98.2 


1.7 / 96.6 
1.7 / 96.6 
1.5 / 95.8 
1.5 / 96.2 
1.7 / 96.6 


VW (original) 
VW (mask-filled) 
Q (mask-filled) 
V (mask-filled) 
W (mask-filled) 


2.3 / 99.7 
2.3 / 99.7 
2.2 / 99.6 
2.2 / 99.6 
2.2 / 99.6 


2.7 / 99.8 
2.3 / 98.7 
2.6 / 99.8 
2.6 / 99.8 
2.6 / 99.8 


1.7 / 97.1 
1.7 / 97.6 

1.6 / 96.6 

1.7 / 97.2 

1.8 / 98.2 


v2 : 
















VW (original) 
VW (mask-filled) 
Q (mask-filled) 
V (mask-filled) 
W (mask-filled) 


2.0 / 95.6 
1.6 / 93.4 
0.7 / 83.1 
1.2 / 90.5 
1.4 / 92.3 


5.5 / 99.7 
4.3 / 99.3 
2.3 / 96.0 
4.0 / 98.9 
2.8 / 96.7 


0.1 / 68.5 
0.3 / 72.9 
1.2 / 89.4 
0.6 / 79.2 
1.9 / 71.6 


VW (original) 
VW (mask-filled) 
Q (mask-filled) 
V (mask-filled) 
W (mask-filled) 


1.3 / 90.7 
1.1 / 85.6 
0.4 / 64.2 
1.1 / 84.9 
1.0 / 83.9 


2.7 / 99.8 

4.3 / 99.3 
1.7 / 95.6 

2.4 / 99.0 
1.9 / 96.4 


0.5 / 68.9 
0.6 / 69.8 
1.1 / 85.8 
0.6 / 71.4 
0.2 / 55.9 



VW (original) 
VW (mask-filled) 
Q (mask-filled) 
V (mask-filled) 
W (mask-filled) 



2.3 / 97.4 
2.1 / 97.1 

1.8 / 96.3 

1.9 / 96.6 
2.0 / 96.4 



4.2 / 99.1 
3.7 / 98.8 
2.9 / 98.3 
3.5 / 98.8 
3.0 / 98.5 



1.3 / 93.5 
1.3 / 94.2 
1.6 / 95.5 
1.3 / 94.1 
1.3 / 93.7 



Table 2. The cr-normalised deviations S and the empirical prob- 
abilities p of the scale-independent diagonal x^-statistics from the 
equations (|6} to ([s} for the different bands and methods as well 
as for the Full Sky and the single hemispheres. 



trum (Hansen, Banday & Gorski 2004b ; Hansen et al. 2008) 



and higher order correlation functions ( ^Eriksen et al. 2004a ) , 
spherical wavelets ( Vielva et al. [2004 ) , local curvature anal- 
ysis (Hansen et al. 2004a), two-dimensional genus measure- 



ments (Park 2004) as well as all three Minkowski func- 



tional ( Eriksen et al.||2004b ), correlated component analy- 



Bonaldi et al. ||2007 ), spherical needlets (Pietrobon et 



2008 



al. 2008), frequentist analysis of the bispectrum (Land & 
Magueijo 2005), two-point correlation functions (Bernui 
Bernui & Hipolito-Ricaldi |[2008) and Bayesian anal- 



ysis of the dipole modulated signal model ( [Hoftuft et al. 



2009). To take a closer look at asymmetries in the WMAP 
five-year data in our investigations, we perform an analy- 
sis of rotated hemispheres, as it was done for the three- 
year data in Rath et al. (2007): For 3072 different angles. 



we rotate the original and simulated maps and then com- 
pute S{r) for the above statistics (mean, standard deviation 
and x^-combination) by only using the data in the result- 
ing new upper hemisphere. Thus, the colour of each pixel 
in the corresponding figure [7] expresses the positive or nega- 
tive cr-normalised deviation S{r) of the hemisphere around 
that pixel in the WMAP-data compared to the hemispheres 
around that pixel in the simulations. We apply this anal- 
ysis for the co-added VW-band as well as for the single 
bands, whereas for the VW-band we use both the original 
and the mask filling method, but for the single bands the 
filling method only. In all charts of figure [7| we can detect an 
obvious asymmetry in the data: The largest deviations be- 
tween the data and the simulations are exclusively obtained 
for rotations pointing to northern directions relative to the 
galactic coordinate system. The maximum value for S{r) of 
the analysis (right column of figure |7| using the mask- 



(0.2)- 



VW (original) 
VW (mask-filled) 
Q (mask-filled) 
V (mask-filled) 
W (mask-filled) 



1.8 / 95.3 
1.7 / 94.5 
1.2 / 90.9 
1.6 / 94.2 
1.5 / 93.7 



5.2 / 99.4 

4.4 / 99.1 

3.4 / 98.6 

4.2 / 99.2 

3.3 / 98.6 



0.5 / 81.6 
0.6 / 82.6 
0.8 / 97.1 
0.5 / 81.9 
0.5 / 81.0 



Table 3. Same as table [2] but for the scale-dependent statistics 
from the equations (|3| to ([sj for the single scale r = 0.2. 



filling method on the co-added VW-band is obtained in the 
reference frame pointing to {0, (j)) — (27°, 35°), which is close 
to the galactic north pole. This proximity to the pole is con- 



sistent to the results of Hansen et al. (2004a) and Rath 



et al. (2007), as well as to those findings of Hansen et al. 
(2004b) and Eriksen et al. ( 2004a) that considers large an- 
gular scales. For the standard deviation (central column of 
figure [7| , the northern and southern hemispheres offer dif- 
ferent algebraic signs. The negative S{r) of the north implies 
a lower variability than the simulations in this region, while 
the south shows a converse behaviour. The fact that the 
plots using the new method show slightly lower values for 
S{r) than the ones using the old method may be explained 
by the fraction of pure noise values within every rotated 
hemisphere, that diminish the degree of difference between 
the data and the simulations. 

Another remarkable feature of figure [T] is the high cor- 
relation between the different bands, that is visible to the 
naked eye but also confirmed mathematically: By calculat- 
ing correlations c among all combinations of those bands 
where the mask filling method was applied, we obtain for 
the mean c ^ 0.99 and for the standard deviation as well 
as the x^-statistics c ^ 0.95. While the Q-band is heav- 
ily foreground- affected, first of all by synchroton radiation 
as well as radiation from electron- ion scattering ("free- free 
emission"), the W-band is mainly distorted by Dust emis- 
sion. The V-band is affected by all these foregrounds, even 
though less than the other bands. As mentioned in chap- 
ter [2] we use the foreground-reduced maps in our analysis, 
but one could still expect some small interferences. Despite 
the different influences on the different bands, we obtain the 
same signatures of non-Gaussianity in all single bands as 
well as in the co-added VW-band. Therefore we conclude 
that the measured asymmetry is very unlikely the result of 
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Figure 7. The cr-normalised deviations S{r) of the rotated hemispheres at the scale parameter r = 0.2 for the mean (left column), the 
standard deviation (central column) and the diagonal x^-statistics (right column) for the co-added VW-Band without (top row) and 
with (second row) the appliance of the mask filling method, as well as for the single Q-, V- and W-bands (third to fifth row), for which 
the mask filling method was always applied. Notice the different colour scaling for each plot. 



a foreground influence but has to be concluded of thermal dices (|Rath et al. 2007) or the Kolmogorov stochasticity 



origin. 



4.2 Local features 



An interesting anomaly in the CMB data is that there 
are small regions which show very high or very low val- 
ues in some local structure analysis. Vielva et al. (2004) de- 
tected the first of these regions, the well-known cold spot at 
(^, (j)) = (147°, 209°) a few years ago by using a wavelet anal- 
ysis. This Spot was re-detected several times using amongst 
others wavelet analysis ( Mukherjee Wang ||2004 Cay on, 
Jin Treas'ter|[2005l |Cruz et al. ||2005| |2006| 7^aling in- 



parameter (Gurzadyan & Kocharyan 2008|. Furthermore, 
there have been some investigations which, in addition to 
the re-detection of the first spot, detected secondary spots 



via directional ([McEwen et al.||2005[ |2006| |2008|) or steer- 
able wavelets ( j Vielva et al. 2007|), needlets (Pietrobon et 



al. |2008|) and again the Kolmogorov stochasticity param 



eter (Gurzadyan et al. 2009). These spots could be the 



result of some yet not fully understood physical process. 
For the cold spot lots of theories already exist which try 
to explain its origin by second-order gravitational effects 
(Tomita 2005 Tomita & Inoue 2008), a finite universe 



model ( Adler, Bjorken Overduin ||2006J , large dust-filled 
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Figure 8. The pixel-wise deviations Si^r of the primal (upper left) and of the smoothed scaling indices map (upper right), both based 
on the VW-band and the scale parameter r = 0.2. The plot in the lower left only shows the values ^ —3.0 of the smoothed method. 
Except for the very small spots in the right part of this mapping, these regions are added to the KQ75-mask. The result is illustrated in 
the lower right plot. 



voids ( llnoue k Silk 12 006", "2007"; 'Ru dnick, Brown k Williams 



2007 



al. 



Granett, Neyrinck Szapudi [2008 ), cosmic textures 
[Cruz et al. | |2007t , non-Gaussian modulation (iNaselski et 



][2008 ) 



2007), topological defects (Battye, Garbrecht & Pilaft- 
20081, textures in a brane world model (Cembranos et 



al. 12008 ) or an asymptotically fiat Lemaitre-Tolman-Bondi 
model ( [Garcia-Bellido Haugbolle |2008| [Masina Notari 



For our investigations concerning spots in the WMAP 
data we only use the mask-filling method of chapter |3.2| 
due to the reasons already explained above. We extend 
the analysis of scaling indices by applying two different ap- 
proaches to detect anomalies: The first one is to calculate 
the cr-normalised deviation of every pixel on the a-response 
of the CMB map. For a given scale parameter r, this is 
achieved by comparing the scaling index a(pi,r) of each 



vector pi^i = 1 , 



of the original data with the mean 



of the corresponding values ai{pi,r),£ — 1, ...,Nsim, of the 
simulations depending on their standard deviation, where 
Nsim denotes the number of the simulations. Formally, this 
reads as: 



(10) 



with 



2 _ 



N, 



N ■ 
Sim 

——- {ai{p^, 

m, J- „ . 



The results are illustrated in the upper left part of figure [S] 
The second approach smoothes the a-maps of the orig- 
inal and simulated data by computing for every pixel the 



mean value of its surroundings given by some specified max- 
imum distance, which equals 3° in our analysis. We apply 
the pixel- wise deviations Si^r again on the resulting maps. 
The outcome of this procedure is shown in the upper right 
part of figure [S] In the lower left plot of the same figure 
only the deviations Si,r ^ —3.0 are illustrated to gain yet 
another clearer view on the interesting areas. 

The first approach clearly shows the cold spot and in- 
dicates some secondary spots in the southern as well as in 
the northern hemisphere. These get confirmed in the plot 
of the smoothing method, where we obtain a deviation of 
up to —7(7 for several clearly visible areas: In the southern 
hemisphere we detect a cold spot at (^,0) = (124°, 320°) 
and another one at (^,0) = (124°, 78°). Both were already 
detected with the above mentioned directional and steerable 
wavelet as well as with a needlet analysis. The former one is 
a hot spot in these investigations. In our analysis, the latter 
spot actually appears as two spots close to each other, which 
is in agreement with jPietrobon et al. | (|2b08 ) . We discover 
another southern cold spot at (0, cj)) — (120°, 155°) which is 
very close to the mask. This spot represents a good example 
for the use of the mask filling method since it is situated at 
the edge of the non-masked region: The influence of the mask 
is diminishing the results of the calculation of the scaling in- 
dices in the area of this spot. This becomes obvious if one re- 
calls the lower left plot of figure^ in which the coordinates of 
the spot would be completely located in a "blue" region with 
low a- values. Since the results of the scaling indices of lo- 
cal features show a similar, namely lower- valued, behaviour, 
an overlapping like that could prevent the detection of such 
spots close to the mask. By using the mask filling method, 
the detection of this cold spot on the edge of the mask is 
equivalent to a detection in an unmasked region, and there- 
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Full Sky 



Full Sky 





Northern Hemisphere 



Northern Hemisphere 





Southern Hemisphere 



Southern Hemisphere 





Figure 9. The cr-normalised deviations of the mask-fining 
method for the original KQ75-mask (blue) and for the modified 
mask of the previous figure (green) in absolute values, plotted 
as a function of the scale parameter, whereby as above "+" de- 
notes the mean, "*" the standard deviation and the boxes the 
X^-combination. The full sky as well as again the single hemi- 
spheres are considered. The blue lines exactly correspond to the 
blue lines of figure [s] 



fore reliable. The spot at (^, 0) = (136°, 173°), described by 
McEwen et ll] 12005) and [Pietrobon et al. | ( |2008| ), is not 
traced in our analysis. In the northern hemisphere, our in- 
vestigation shows two other cold spots at (^, 0) = (49°, 245°) 
and (^,0) — (68°, 204°), which do not correspond with the 



so-called northern cold spot of Gurzadyan et al. (2009), but 



with the results of McEwen et al. (2005), where again one 



of them is a hot spot. Also Pietrobon et al. (2008) locates 



one of these two spots. All these results were achieved with 
an analysis of the VW-band, but we find similar results in 
a single band analysis. 

It is possible to define a new coordinate frame, including 
a new direction of the "north pole", such that all of these 
spots are contained in the "southern" hemisphere. This new 
north pole would then be located at (^, 0) = (51°, 21°). 

If the considered spots really depend on some yet not 
completely understood, maybe secondary, physical effect, 
they should not be implemented in a testing for intrinsic 
non-Gaussianity. For this reason, we modifiy the KQ75-mask 
by additionally excluding all above mentioned spots. A small 
peculiarity at the edge of the mask next to the cold spot as 



Figure 10. Same as figure [ojbut exclusive for the modified mask 
of figure [s] applied to the Q- (red), V- (green) and W-band (blue). 
This plot is associated with figure |6] 



well as three very small blurs in the right half of the lower left 
mollweide projection in figure [S] are not considered, since we 
regard their appearence as insufficent for being a distinctive 
feature. The modification of the KQ75-mask is illustrated in 
the lower right part of figure [S] 

We now apply this new mask to the a-response of both 
the WMAP data as well as the simulations and repeat the 
analysis of chapter [4. 1| The results are illustrated in the fig- 
ures [9] and 10 as well as in table [4] A clear increase of S{r) in 
comparison to the former analysis is evident. This height- 
ening is in particular present in the southern hemisphere, 
where we detected more local features than in the north. The 
largest increase takes place in the co-added VW-band, where 
we now reach deviations of up to 4.0 for the x^-combination 
in a full-sky analysis (former maximum: 2.9) and to the ex- 
tend of 6.0 in an analysis of the northern hemisphere (for- 
mer maximum: 5.5). But also the single bands in figure 10 
as well as all scale-independent diagonal -statistics in ta- 
ble |4] show without exception a greater evidence for non- 
Gaussianity. 

One could have expected to obtain higher values for 
S{r) since the a-response of the WMAP data in comparison 
to the one of the simulations featured a shift to higher values 
(see figure [3|: By now cutting out the local features, that 
exclusively consist of cold spots in terms of pixel- wise devia- 
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Full Sky 


Northern Sky 


Southern Sky 




(S/yo) 






VW fmask-fillprl'l 


2.4 / 97.6 


2 8 / 98 4 


2 / 97 2 


Q (mask-filled) 


2.3 / 97.5 


2.8 / 98.3 


1.8 / 96.9 


V (mask-filled) 


2.3 / 97.6 


2.7 / 98.1 


1.9 / 97.1 


W (mask-filled) 


2.4 / 97.7 


2.6 / 98.2 


2.0 / 97.4 


VW (mask-filled) 


2.6 / 96.7 


4.8 / 99.8 


0.2 / 97.2 


Q (mask-filled) 


1.3 / 90.6 


2.5 / 97.0 


0.6 / 82.1 


V (mask-filled) 


2.0 / 94.9 


4.4 / 99.4 


0.4 / 74.1 


W (mask-filled) 


2.2 / 96.2 


3.1 / 98.0 


0.4 / 78.9 


VW (mask-filled) 


2.7 / 98.0 


4.0 / 99.1 


1.6 / 96.2 


Q (mask-filled) 


2.2 / 97.1 


3.0 / 98.6 


1.6 / 95.6 


V (mask-filled) 


2.4 / 97.3 


3.7 / 98.9 


1.5 / 95.3 


W (mask-filled) 


2.5 / 97.9 


3.2 / 98.7 


1.6 / 95.5 


Table 4. Same as table [2] but after excluding the cold spots via 
the modified KQ75-mask. 

tions, one excludes spots that showed lower values than the 



average of the simulations (see equation 10). Therefore, the 
shift to higher values becomes even larger, hence leading to 
a higher S{r). Still, the exclusion of the spots is helpful and 
necessary, since these local anomalies could origin in some 
independend physical process, as mentioned above. 



5 SUMMARY 

We performed a scaling index analysis of the WMAP 5-year 



data following up the investigations of Rath et al. (2007). 
For more realistic results around the mask, we additionally 
implemented a mask- filling method. By comparing the Q-, 
V-, W- and the co- added VW-band of the WMAP data with 
1000 simulated maps per band, we (re) detected strong devi- 
ations from Gaussianity as well as asymmetries in the data, 
which can be summarized and interpreted as follows: 

The scaling index values of the WMAP data are shifted 
to higher values and feature a higher variability than those of 
the simulations, especially in the northern hemisphere. This 
effect can be interpreted as less structure as well as more 
structural variations in the CMB signal compared to the 
corresponding Gaussian model. The results are confirmed 
by several statistics, that show deviations from Gaussianity 
of up to 5.9(7 in the scale-dependent, and 5.5cr in the scale- 
independent case. These results are slightly lower applying 
the mask- filling method, and show high similarities within 
the different bands. In addition, we detected strong asym- 
metries by performing an analysis of rotated hemispheres: 
rotations pointing to northern directions show by far higher 
higher deviations from Gaussianity for the mean and the 
X^-analysis than rotations pointing to the south. Observing 
the standard deviation, we obtained a negative outcome in 
the north and a positive in the south. This implies that the 
north possesses a more consistent pattern than the simula- 
tions, while the south shows the converse behaviour. This 



feature is in line with later investigations of local features, 
where we detected more local anomalies in the southern than 
in the northern hemisphere. 

Furthermore, we performed an analysis of local features 
by studying pixel-wise deviations from Gaussianity with 
and without a previous smoothing of the a-responses. For 
these investigations, we exclusively applied the mask filling 
method which can reduce the distorting effects on measure- 
ments like the scaling indices that appear when cutting out 
the masked regions. This mask-filling method eliminates the 
diluting effects on the border and therefore allows for an 
analysis of local features, which show a similar behaviour of 
lower outcome in the scaling index method. We detected the 
well-known cold spot and three additional spots in the south- 
ern as well as two spots in the northern hemisphere. Except 
for one single spot in the south, all findings are in agree- 
ment with former results of different investigations. Since 
these spots could origin on some yet not completely under- 
stood physical effect, we excluded them from the data set 
and repeated the former analysis. Instead of obtaining lower 
deviations, the results show an increase of no n- Gaussianity 
in all bands. Therefore, the discovered local anomalies are 
not the reason of the global detection of no n- Gaussianity, 
but were actually dampening the deviations on average. In 
former isotropic wavelet-based analyses, an exclusion of de- 
tected spots lessened the significance level of indications of 
non- Gaussianity (Vielva et al. | [20Q4). Our new findings in- 
dicate in contrast, that the isotropic scaling index method 
can detect several different yet complementary aspects of 
the structural composition of the underlying data. The re- 
sults of our investigation are in agreement with the steer- 
able wavelet-based analysis in Wiaux et al. ( 2008| ), where 
the non-Gaussianites were conserved after excluding the de- 
tected local anomalies. 



6 CONCLUSIONS 

The redetection of indications for non- Gaussianity of the 
WMAP 3-year data analysis leads to the conclution that 
the observed results are not time-depending. In contrary, 
we can detect even higher deviations from the simulations 
which mimic the Gaussian properties of the best fit ACDM- 
model. Therefore, it is highly improbable for the results to 
be caused by effects related to short-term measurements. 

In addition, the coherence between the different anal- 
ysed bands implies that the foreground influence plays only 
a minor role but that the results are very unlikely to be 
truely of thermal origin. 

Finally, the agreement of the detected spots with former 
investigations confirms the existence of these local anoma- 
lies. 

The two most important tasks for future studies are: 
First, to identify possible reasons for the indications of non- 
Gaussianity, which could be possible with the attainment 
of more and more precise data, e.g. with the upcoming 
PLANCK-mission. Second, to figure out possible sources of 
the observed local features and thereby solving the question, 
if these anomalies are due to systematics or foreground ef- 
fects or indeed represent variations in the CMB signal itself. 
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